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SPIN EFFECTS IN GRAVITATIONAL RADIATION BACKREACTION 
III. COMPACT BINARIES WITH TWO SPINNING COMPONENTS 
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KFKI Research Institute for Particle and Nuclear 
Physics, Budapest II 4 , P.O.Box 49, H-1525 Hungary 


The secular evolution of a spinning, massive binary system in eccentric orbit is analyzed, expand¬ 
ing and generalizing our previous treatments of the Lense-Thirring motion and the one-spin limit. 
The spin-orbit and spin-spin effects up to the 3/2 post-Newtonian order are considered, both in the 
equations of motion and in the radiative losses. The description of the orbit in terms of the true 
anomaly parametrization provides a simple averaging technique, based on the residue theorem, over 
eccentric orbits. The evolution equations of the angle variables characterizing the relative orienta¬ 
tion of the spin and orbital angular momenta reveal a speed-up effect due to the eccentricity. The 
dissipative evolution of the relevant dynamical and angular variables is presented in the form of a 
closed system of differential equations. 


I. INTRODUCTION 

The motion of a binary system of two spinning bodies under the influence of gravitational radiation has been treated 
in several recent works. The aim of these investigations is to obtain a self-consistent description of the evolution, such 
that templates for the gravitational radiation pattern can be provided for the wave observatories under construction. 
Spin effects modulate both the amplitude and the frequency of the waves. As of now, the basic properties of the 
radiative evolution have been understood in a perturbative framework, both in a post-Newtonian expansion |l]- and 
by black-hole perturbation techniques i- i- The post-Newtonian expansion proceeds in powers of the parameter 
e w w Gm/c^r. In both approaches, the radiative losses in characteristic quantities, with the inclusion of 

spin-orbit and spin-spin effects, have been computed. The averages over circular orbits have been obtained in Q] and 

m in 

The importance of eccentric orbits in various physical scenarios has been emphasized by several authors [1^]- [[13| . 
Quinlan and Shapiro 0 argue that clusters ingalactic nuclei in the final stage of collapse will contain a significant 
number of eccentric binaries. Hills and Bender suggest that many massive (M « 10® — lO^M©) compact objects in 
the galactic centers are gravitationally deflected by others and there is insufficient time left for circularization before 
plunging. The behavior of eccentric binaries under radiation reaction forces has been studied in 0 , i and 0 , and 
(neglecting spin effects) in |0. Averaged radiative losses for eccentric orbits have been obtained for a test particle 
by Ryan 0. For finite masses, partial descriptions (by computing the radiation losses of E and L) have been given 
by Rieth and Schafer |^. The averaging procedure for eccentric orbits has proven cumbersome in these works. 

In the present series of papers, we investigate the influence of the intrinsic spin on the evolution of a radiating 
eccentric binary system in a post-Newtonian approach. The presence of spins complicates the description of the 
orbit considerabN; indeed the computations neglecting spins have reached precisions up to the 7/2 PN order |p^ - 
0. In papers 0 and (to be referred to as I and II, respectively) we have chiseled two convenient tools for our 
investigation. The first is the parametrization of the orbit by the generalized eccentric anomaly (^) and the true 
anomaly (y) parameters. These parametrizations have already been employed in various special cases. The test- 
particle limit has been considered in I and the one-spin limit in II. The second ingredient of our technology is the use 
of the residue theorem for averaging the gravitational radiation losses by means of these parametrizations. We now 
further develop and employ both of these tools for the treatment of two spinning masses. 

For the purpose of predicting the evolution of a binary with two spins, it is crucial to determine the variation of 
a complete set of geometrical parameters characterizing not just the orbit but also the orientation of the angular 
momenta. The amplitude and polarization angle of the gravitational signal is modulated (as has been stressed in 0 ) 
by the changing orientation of the binary in the observer’s frame of reference. 

In this paper we characterize the motion by computing the evolution of the system parameters to the 3/2 PN order, 
including spin-orbit and (as an order-of-magnitude analysis will reveal) spin-spin contributions. 

The plan of the paper is as follows. In Sec. II., the motion of the binary system is described following |^, 0 and 
1^. We use the covariant spin supplementary condition {SSC) 1^. Then we introduce the conserved quantities of the 


I 


nonradiative motion: the energy E and the magnitude of orbital angular momentum L, further the angles subtended 
by the orbital angular momenta and spins. The evolution of these angles has been discussed for circular orbits in [^ , 
[po| and Q. Here we give both the instantaneous and averaged evolution equations of these angles, to 3/2 PN order, 
including the relevant spin-spin terms, for eccentric orbits. 

In Sec. Ill we compute the instantaneous losses of the energy i?, magnitude L of the orbital angular momentum 


and the angle variables from Kidder’s universal expressions (in [g), and by use of the Burke-Thorne ||^ potential. 

We derive the secular evolution of the system variables due to gravitational radiation reaction (Sec. IV). The 
contributions from the Burke-Thorne potential have no effect on the averages over a period of revolution. Our results 
provide a dual desc ri ptio n [Eqs. (EII )-€1 in terms of the energy E, orbital angular momentum L and the spin 
angles vs. Eqs. (4.S)- (|4.9| ) in terms of the semimajor axis a, eccentricity e and the spin angles] of the radiating binary 
system. The results in I and II are special limiting cases. 


II. THE ORBIT OF THE BINARY SYSTEM 


The bound state of a two-body system with masses mi and m 2 and spins Si and S 2 , respectively, is described by 
the Lagrangian 


C = 


Gm^i 2G/r 


[r X (S -I- cr)] -I- 




(a X cr) . 


( 2 . 1 ) 


2 ' r ' '■ \ ' /j ' 2c^m 

Here r = jrj is the relative distance, v the relative velocity, fi the reduced mass and m the total mass of the system. 


77117712 7712 

fl = - , 777 = 7771 -I- 7772 , V = - ■ 

mi 


( 2 . 2 ) 


TTll -I- 7772 

The total and weighted spins are defined 

S = Si-bS 2 , cr = 7781 + 77-^82 . 

We have kept only the leading-order spin-orbit coupling, and adopted the spin supplementary condition of The 
relative acceleration entering the Lagrangian is 


(2.3) 


Gm G ( r 

a = --I- -^r [(r X v) • (8 -|- cr)] - v x (48 + 3cr) -|- 3-r x (28 -|- cr) 


and the momenta: 


dC 77 

Q = = —7;—cr X V 

dsL 2c^m 


(2.4) 


(2.5) 


dC . G/7 

P = ^-q = 7iv-h X (28 + cr) . 

av c^r'^ 


( 2 . 6 ) 


Here an overdot denotes derivative with respect to the time parameter t. The constants of motion are the energy 
E = p -v + q a. — C and the total angular momentum J = 8 -|- L. The orbital angular momentum is 


L = rxp-|-vxq = LN + Lso , 

with the Newtonian orbital angular momentum and spin-orbit terms 

Ln = /ir X V , 


Lso = 


^ ^ + <^))] - ^ [v X (v X ,t)]| . 


(2.7) 

( 2 . 8 ) 

(2.9) 


The spin precession equations have been given in 

G /A + St] 


Si = 


LN-S2 + ^(r-82)r x 81 , 


S 2 = ^ f ' ln - 81 + ^(r . 8 i)r ) x 82 . 


( 2 . 10 ) 
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The magnitudes Si and S 2 of both spins are separately constants. 

All spin effects are formally of 1 PN order. A simple estimate S/L « reveals however that the spin -orbi t 
and spin-spin terms are of respective orders and in general. Exceptionally, in the spin precession Eqs. (2.1C) 
(predicting the evolution of the directions Sj of spins), the spin-orbit and spin-spin contributions are of order e and 

g3/2_ 

The precession rate of the orbital angular momentum, however, will not inclu de spin-spin terms, as they are of 
order. From the conservation of the total angular momentum J and from (2.10| ) we obtain 


L = -S = 


G 


2 c^r^ 


(4S -C 3cr) X L . 


( 2 . 11 ) 


We substituted here, as we did in II, L in place of the Newtonian angular momentum Ln. It follows from ( 2.11 ) that 
the magnitude L of the orbital momentum and its projection on the vector 4S -|- 3<t are conserved. However, the total 
spin S shows a more complicated motion pattern. 

The energy and the orbital momentum square are 


Gm^i ^ G'(L-(t) , ^ 2/^2 , „;„ 2 a,;, 2 M Gmn ^ G(L-cr) 


2 r 

= /rV^(02 -Csin^e (^2) -4 


, 2,3 =^[f^ + r\0^+sin^9^^)]- 


G/x(L • S) 2 A(L-cr) 


( 2 . 12 ) 

(2.13) 


Hence we express and in terms of r, constants of the motion, and in terms of L-cr and L • S as follows. 


v=-E- 


2Gm 2G(L-cr) 




• 2 r.-® 

’’’ ~ 2 -^ 2 -TT 

fj, r 


-f 2 


A(L-cr) 2G 




(2L • S 4 -L-cr) . 


(2.14) 

(2.15) 


Even though the spin projections of L are not conserved, only the leading-order (conserved) contributions from the 
quantities L • Si appear in Eqs. (2.14) and ( 2.15| ). 

To complete the description of the motion, we next compute the evolution of the angle 7 subtended by the spins 
Si and S 2 , further the angles Ki subtended by Si and L. In terms of these angles we may express 


L • S = L{Si cos Ki -I- S 2 cos K 2 ) , L'cr = L{r]Si cos ki + rj ^82 cos K 2 ) ■ 

The variation of the angles Ki and 7 is obtained from Eqs. o of the spin precession: 

, 3G{2 + r]-^) 82 ^ A ^ 3G 52, A 

(cosKi) =---L • (Si X S2)+^-7- —(r • S2)L • (r x Si) , 


(2.16) 


2^2 J ,3 


p2j,5 


, ,. 3G(2 -|- rj) Si T An 3G 5i . £> \t i £> \ 

(COSK2) =--pL • (Si X S2)+^-7- —(r • Si)L • (r X S2) 


2c^r^ 


L 


(2.17) 


, 3G(7-7?-1)^ ~ , 3G(r-S2-r-Si) - . 

(COS 7 ) =-— -E ■ (Si X S 2 )-f - 2 ”^- -Y ■ (Si X S 2 ) . 


All terms are of order excepting the first term in (cos 7 )', which is of order e. Thus in Eqs. ( ^.17| ) we need the 
mixed products only to leading order. First we evaluate 


L- (Si X S2) =p 


(r-Si)(vS2)-(r-S2)(v-Si) 


L • (r X Si) 


/ir 


r(v • Si) -r(r • Si) 


(2.18) 


The scalar products may be expressed by use of the angles V' and ipi, subtended by the node line (cf. Fig.l in II) with 
the momentary position vector and the respective projections of the spins in the plane of orbit. These expressions 
are: 


r • Si = r sin Ki cos('(/' — ipi) , 

V • Si = r sin Ki cos{ip — ipi) -sin Ki sm{ifj — ijji) . 


(2.19) 
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We have used here the time derivative of the angle i.e., ip = L/ 

The angles ipi do not independently vary. This can be seen by noting that the node line is orthogonal both to the 
total angular momentum J and to the Newtonian orbital momentum Ln- The projection of the total spin on the 
node line yields 


Si sin m cos ipi + S 2 sin K 2 cos ip 2 = 0 ■ ( 2 . 20 ) 

Hence we may express both angles ipi in terms of, say, l^ip = ip 2 — ipi- Furthermore, from the spherical cosine identity 
the angle Aip is determined by the angles Ki and 7 (Fig.l): 

cos 7 = cos Ki cos K2 + COS Alp sin ki sin K2 ■ (2-21) 



So 


FIG. 1. The angles subtended by Ln and Si , and the difference Aip of the angles subtended by the node line and the spin 
projections on the plane of orbit. To leading order, the Newtonian angular momentum Ln = L. These angles are related by 
the spherical triangle identity, ip and ipo are the respective angles of the position r and the direction of the periastron with the 
node line. 


Thus, in terms of the angles Ki and Aip the expressions (2.18), to leading order, take the form: 

L • (Si X S 2 ) = L sin ki sin K 2 sin Aip, 

L • (r X Sj) = —Lr sin Ki sin {ip — Ipi). 

The last mixed product in ( ^.17 ), again to leading order, works out as 

r • (Si X S2) = r[cos ki sinK2 sin{ip — 1 P 2 ) — cos K 2 sin/ti sin{ip — ipi)] . 


( 2 . 22 ) 

(2.23) 


It should be stressed that our treatment of the angles is fully coordinate invariant. 

We proceed next with the parametrization of the orbit. Following the approach described in Sec. Ill of II, the 
eccentric anomaly parameter ^ is introduced by 


Gmpi Gfr(2L • S + L-cr ) 

2 F; ^ 


^0 

2E 


G2to^ 2(2L. S + L-cr) E(L-a-) 




0 


c'^mp.Ao 


cos 5 


where Aq is the length of the Runge-Lenz vector to the zeroth order in the spin: 

- 2 \ 1/2 


Ho = + 


2EL^ 

As shown in II, the orbital period is given by the integral of dt/d^ from 0 to 27r, 

GmpP 


T = 2tt- 


{-2fiE)i 


(2.24) 


(2.25) 


(2.26) 


4 

















The true anomaly parameter x with properties explained in I (and II) is introduced similarly, 

2 G(2L-S + L-cr) Ao(2G^m'^+ EL^) + Gmfi(2G^m^+ 3EL^) cosx 

T =■ - -|- - - 

p,{GmpL + Aq cos x) c^L'^Ao {Gmfi + Aq cos x)^ 

2E{L-cr) Gmp^AQ + {G‘^m?p^ + EL'^) cos x 


c^mfEAo 


{Gmp, + Aq cosx)^ 


(2.27) 


This yields 


dt 1 dj' LL'T^ r 1 1 

— = T-;— = -r- “I 1-^ 5 — T 7 [(SL • S + E-(T)Gmp^iZGmp + Aq cosy) — EL^{L-(t)\ > . (2.28a) 

ax r ax E [ c^mL^ ■' J 


The leading order expressions of the angle variable tp and of r are also needed in parameterized form: 

= V'o + X , 


^0 . 

r = -r- sin X 


L 


(2.29) 


With ( ^.27 ) and (2.29), all expressions we are interested in are expressed in terms of the true anomaly parameter x- 
We integrate these expressions as follows. 


/»27r Ti 

E{t)dt = / Fix)— dx . 

Jo dx 


(2.30) 


The integration is carried out by computing the residues enclosed in the circle ^ = e*^. As a rule, we find that there 
is only one pole0, at = 0. 

For a first application, we compute the average rate of change of the angles. 


,, 3Gil + p-^)i-2 pE)^/^S2 . . . , , 

((cosACi) ) =--sin Ati sin a ^2 sin , 

2 c^L'^ 

,, 3Gi\Ev0i-2pEfl'^S^ . . . , , 

((cos K 2 )) = --sin m sm /^2 sin iX'ijj , 

Zc^L’^ 


(2.31) 


((cosy)’ ) = 


3Gi-2pE)^/'^ 

2 c 2 L 2 


^2 , . . ... 

77 — 77 + — cos Ki -— COS K2 sm Ki sm K 2 sm i\ip . 


L 


L 


The secular changes do not vanish. This should be contrasted with the behavior of the one-spin system, presented in 
II, where the angular momenta are frozen in a rigidly rotating parallelogram, and the relative angles are constant. 

We may express the secular changes of the angles in terms of the semimajor axis a = {Gmp) /(—2i?) and eccentricity 
e given by 1 — = {—2EL'^)/{G^mJp^). The Keplerian values are sufficient to the accuracy needed here. We then 

get 


((cOSKi)’) 
((cos K2)') 

((cosy)’ ) = 


3G(1 + 77-1)^2 
2c2a^(l — e2)3/2 
3 G (1 +77)^1 


2c‘^a^(l — e2)3/2 


sin Ki sin K 2 sin Atp , 
sin Ki sin K 2 sin Atp 


(2.32) 


3 G 


1 /2 

(77 — r]~^)p (G77ia(l — e^)) -I- Si coski — S2 COSK2 


2c2a^(l — e2)3/2 


sin Ki sin K 2 sin Atp . 


In the particular case when the masses and spin magnitudes are equal, mi = m 2 and Si = S 2 and with the 
appropriate change of notation, these equations agree with Apostolatos’ Q Eqs. (4) up to a constant multiplier 
1/(1 — of the time t. Thus the eccentricity e of the orbit speeds up the evolution of the angular variables Ki 

and y. In the generic case, with arbitrary values of the masses and spin magnitudes, the scaling factor is still present 
in the denominators of Eqs. ( p.31 ), but no corresponding term in the numerators. We thus find that the eccentricity 
of the orbit accelerates the evolution of the spin directions. 


We were initially unaware of the conditions under which the true anomaly parameter has this property. A systematic approach 
has recently been developed to the parametrizations of the perturbed Kepler motion. By its use, we can now prove that 
the poles of these integrands are in the origin of the complex parameter plane. 
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III. INSTANTANEOUS RADIATIVE LOSSES 


We obtain the instantaneous losses of the constants of motion E and L and the angles Ki and 7 using Kiddy’s Q| 
results and the Burke-Thorne potential The variation of the remaining angles follows by the relations (2.2C) and 

(Ell- 

First we find the radiative spin losses. To lowest order the radiation-reaction potential is the Burke-Thorne potential 




(3.1) 


where is the fifth time derivative of the system’s quadrupole-moment tensor and are Cartesian coordinates 
centered on the spinning body. 

The radiative spin loss evaluated in II holds for each of the spins in its system of principal axes of inertia: 


1 d(Si) 
Sj dt 


2G /0 


V 0 ' 


-1 e^.p/(^;(S 0 p(S 0 . 


(3.2) 


Here 0i and 0' are the principal moments of inertia and is the angular velocity of the spinning axisymmetric 
body. These quantities are related by Si = 0'Hi. 

Two important properties of the radiative spin loss dSi/dt emerge from (3.2): 

(z) they are of 2"^^ post-Newtonian order and 

(zz) are each perpendicular to the respective spins Si. 

Next we consider the energy loss dE/dt and total angular momentum loss dJ/dt. These have been computed 
by Kidder using the Blanchet-Damour-Iyer formalism pM. Keeping the Newtonian and spin-orbit terms in the 
expressions, we have 


dE 

8G^m^/z" 


dt 

15c^r‘* 

u 


8G^m/z 

(T 


15c^r® 


dJ 

8 G'^my, 

dt 

5 c°r-^ 



(12^2 - llr^) 

Ln-S) (27f2 -37 z;2 - 12 


(3.3) 


Gm 


+ (Ln-ct) ( 51r^ — 43z;^ -I- 4 


Gm 


Ln -3r^ + 2v^ + 2 


Gm 


4GV^ 


2 Gm 


. Gm 


5 


3 p '(^ “ )(® - X [v X (7S -k 5cr)] 


Gm 


-r X 


(r X S) ( 6r 


17 2 r^C!m\ , N o 2 2 Gm 

— —v^ + 2 - M- (r X (t) 9r^ - 8v^ — - 


r 

--V X 

r 


(r X S) -3072 -k 24v^ + 


r J 
29 Gm 


3 r 


-kv X 

Ln 


2 10 2 23Gm\ 


(v X S) 18r^ - 12z;^ - 


3 r 


r/ ^ I r -2 .2 ^ Gto 

-k 5(r X (t) -5u -k 4z;^ -k -- 

3 r 

35 2 «Gm 


(3.4) 




(Ln • S) 30r^ - 18z;^ - 


3 r / 
92 Gm 


-k (v X cr) ISr'^ - — 9 

3 r 


3 r 


-k (Ln • (t) I 35r^ — 19w^ — 


2 71 Gm\ 


3 r 




The loss dL/dt = h-dJ/dt in the magnitude of the orbital angular momentum follows from (3.4) as the spin-orbit 
terms in dJ/dt do not receive contributions from the radiative spin losses, cf. property (z). All mixed vector products 
in (3.4), when projected to L, can be converted to one of (L ■ S) an d (L • cr). The first, Newtonian, term is expressed 
by use of Eqs. (2.7) and (pT|). Using then Eqs. (2.14), (2.15) and ( pl 6 | ) we obtain purely radial expressions for the 
instantaneous losses: 


dE 

dt 


15c®r® 


(2/iifr^ -k 2Gmfi^r + IIL^) -k 


8G^mL 

15c^/xr® 


(3.5) 


X (20yEr'^ — 12Gmy?r -k 27L^) (^i cos zti -k S 2 cos K 2 ) 

+ (GfiEr'^ — 18Gm/z^r -k 51L^) (zyS'i cos zti -k ??~^S '2 cos K 2 ) 
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(3.6) 


{ 2 fiEr'^ - 3L^) + 


8 G 2 


dL SG^mL 
dt 5c^iir^ ^ ^ ' 15c7^2^7 


X |/7r [l2Gmfi^Er'^ + 3 + 6 EL^) r — llGm^L^] (S'! cos ki + 52 COSK 2 ) 

+ \2y,'^E'^r‘^ + \2Gm^^Er^ + ‘ipL{G^w?y? + ?)EL^)r^ — bGmy^L^r + 15L^] {rjSi coski + ? 7“^52 cos K 2 )| ■ 

We also need the projections of the orbital angular momentum loss in the directions of the spins for evaluating the 
radiative change in the angle variables. The foregoing considerations yield these projections, after a straightforward 
but cumbersome computation, in the form 


dL ^ dL 

—-Si = — COSKj 
dt dt 


2 G^ 


15c^ 


(2pr 

i=i 


„3-2i 


Qij) SjSi 


sin K,i sin hij 


(3.7) 


with 


Pij = 9yL{—2yEr^ + 5L^)rr sin('!/;j — i/ji) 

+ [12Gmy^Er^ + Sy (G'^m^y^ + ISTIL^) + lAGmy^L'^r — cos{'ipj — ipi) 

—i[AGmy^Er^ + y {G^m^y^ + 6EL^) + Gmy^L^r + cos(2V’ — tpi ~ ’’Pj) 

+3yL{6yEr‘^ + 2Gmy^r — 3L^)rrsin(2'0 — ipi — ipj) , (3.8) 

Qij = lbyL{—2yEr^ + 5L^)rr sm{'tpj — ipi) 

EpiyPE'^r'^ + 2AGmyPEr^ + 6^ {G'^mPy^ + 15£1L^) + 29G'my^Lpr — 75L'*] cos{tjjj — ipi) 

-[Ay^E^r^ + 2AGmy^Er^ + 2y {SG^rn^y^ + 13EL^) + SGmy^L^r + 15L^] cos(2V' --ipi - ipj) 

+yL(3AyEr^ + 12Gmy^r — 15L^)rr sin(2^/! —-ipi — ■tpj). (3.9) 


These expressions depend on r as well as on the angular variable V” and on the derivative r. 
Thus the spin projection of the instantaneous change in L is 


dL 

dt 


•Si 


1 

L 


dL 

dt 


•Si 


dL 

-cos«. 


2 G^ 

Ibc^ y^Lr'^ 


2 

{ 2 pij + Sj sin/Ci sinKj . 

i=i 


(3.10) 


Due to property (ii) the instantaneous change in Si equals the right hand side of (3.2) and its projection to L 
yields: 


L 


dSi 


2 G^r 


0 ,: 


>' 1 I • z 

7 - 1 1 sm Ki 


(3.11) 


dt bd’yPr'^VLi \0 

X \AL{—ld>Eyr^ — 20Gmy^r + 15L^) cos(2'0 — 2'0i) — yrr{V2Eyr'^ + 2QGmy^r + 45L^) sin(2'!/) — 2^/;^)] . 
We are now in position to compute the radiative losses of cos Ki : 


d cos Kj _ d - A X _ c^L A , f dSi 
dt dt ^ dt ^ dt 


(3.12) 


Finally we find from Eq. (^) the radiative change in the angle 7 subtended by the spin vectors: 
dcosy 


m 


dt 




?)C^y?r'^ ^ Lli V 0' 


X |/irr {12Eyr^ + 20Gmy^r + 45L^) (sinK^ cos kj sin(2^ — 2il)i) — cos Ki sin/tj sin(2^/! — -0^ —'4’j)) 


+ {12Eyr^ + 20Gmy^r — 15L^) cos Ki sin/Cj sin(0j — ipi) 


(3.13) 


+AL{18Eyr^ + 20Gmy^r — 15L^) (sinK^ cos Kj cos( 20 — 20^) — cos Ki sin Kj cos(20 — 0^ — 0j))|. 

We emphasize that in spite of the fact that the radiative spin losses dSi/dt are of 2""^ post-Newtonian order, they 
contribute at order to the instantaneous angular losses. 
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IV. AVERAGED RADIATIVE LOSSES 


The instantan eous lo sses o f the constants of motion and the angles subtended by the orbital and spin angular 
mome nta, Eqs. ( 3.5 ) - ( 3.13 ) are in a form suitable for parametrization with the true anomaly parameter y, using 
(2.27) and (2.29). Then the averaged losses are computed by the residue theorem, passing to the complex variable 
C = 

The averaging procedure yields for the constants of motion: 


dE 

dt 


dL 

dt 


X { (520E3l® + l07A0G‘^m^^i^E^L'^ + 2A990G‘^m‘^n^EL^ + 12579G®mV®)(-S'i cos ki + S 2 cos K 2 ) 

+ {256E^L^ + 6660G^m^E^+ 16660G‘^m‘^^^EL^ + 867SG^m^^^){r]Si coski +r]~^S 2 cos K 2 )} (4.1) 

4G2m(-2£;/r)3/2 GH-2Efj,)^/^ 

=-> + 15cV> 

X {(IISSE^L"* + 6756G^m^/x^£'L^ + 5345G'^m'^/x®)(S'i cos ki + S2 cos K2) 

+ {772E^L'^ + AA76G^m‘^fi^EL^ + 8665G'^m‘^fi^){r]Si coski + r]~^S2 cos/C 2 )} 

For the angle variables, we first find that ( ^.11 ) and ( 3.13| ) average out: 

(^ 4 )-. 4 )-- 

Thus, when averaging Eq. ( |3.12|) , only the first term survives, which leaves us with 

'dKi\ _ G^(-2E/r)3 /2 
dt / 30c^L® 

+ (884E^L^ + 326AG^m^^^EL^ + 1795G"^TO^/i®)(?7S'i sinwi + ri~^S 2 sinK 2 cos A?/>) 

+6(52E^L^ + 92G^m^EL"^ + 8SG^m‘^fj,^)[Si sinKi cos(2tjji — 22jjo) + S2 sinK2 cos(^/>i+V'2 —2'0o)] 

+2{238E^L'^+ ASlG^rn^fj,^EL\l56G^m‘^fj,^)[r]Sisin KiCOs(2tjji—24’o)+r]~^S2 sinfi;2 cos(^/>i+^/>2— 2 V'o)]} • 

The corresponding expression {dK 2 /dt) is obtained by swapping the indices 1 2 in Q and substituting ry ^ 77 

For completeness we give the radiative evolution equations also in terms of orbital elements suited for our perturba¬ 
tive treatment. The generalized semimajor axis a and eccentricity e are introduced by r^ax = a(l ± e). (The turning 

points r^a.x follow from (2.24) inserting ^ = tt and 0 respectively.) 


(4.2) 

(4.3) 

{(1140E^L^-|-4164Gm^/r^EL^ -I- 2285G^TO^/i®)(S'i sin ki -|-S '2 sin K 2 cos At/i) (4.4) 


1 - 


1 - = - 


2 EL^ 




Gnifi 
2E L 
4(2L • S-f L-cr) 


2E(2L • SL-cr) 


2 mL 2 




2E(L-cr) 


c^mL^ . 


(4.5) 


The inverse relations are 

Gm^ 


E= - 


2 a 


1 -b 


G^/^ (2-I-?y)S'i cosKi-b (2-I-ry ^)5'2C0SK2 


2 G 1/2 


a3/2(l_g2)l/2 

= Gm/r^a(l - e^)|l - ^2^i/2q3/2(x _ e2)3/2 ['^'i cos «;i(3-b -b 2?y) -b S '2 cos s;2(3-b -b 2ry“^)] | . 


(4.6) 


The first terms on the right- hand sides represent the Keplerian approximation. We emphasize here that the expressions 
of the angular precessions ( ^.32| ) are unchanged if we rewrite them in terms of a and e defined above. 

A straightforward computation yields the averaged radiation losses of a, e and Ki'. 


da\ _ 2G^m^fi{37e^ + 2926^ -b 96) 

dt / 15c®a^(l — 


q7I2 3/2 , 

15c7a9/2(l-e2)5 2128)(S'! coski -b S'acos ^ 2 ) 

(291e® -b 42246^ -b 7924e^ -b 1680)(ryS"! cos ki -b ry"^S '2 coskz)} , 


(4.7) 












































de\ (121e^ + 304) 

dt / 15c®a^(l — 

(^7/2 3/2 „ , 

+ 3Q^7qii/2(']^ _ g2)4 + 5592e^ + 7032)(5'i coski + S 2 COSK 2 ) 

+ (1097e'‘ + 6822e^ + 6200)(r7S'i cos m + r]-^S 2 COSK2)} , 


/^7/2 3/2 

= 3Q^7^ii/2(']^ _ g2)4 {(2856^^ + 1512e^ + 488 ){Si sin Ki + S 2 sin K 2 cos Ai/-) 

+ (221e"^ + 1190e^ + 384)(r7S'i sin/ti + ri~^S 2 sinK2 cos Ai/i) 

+ (156e"^ + 240e^) [S'! sin m cos(2//!i — 2'ipii) + S 2 sin k ,2 cos{'ipi + '02 — 20o)] 

+ (119e"^ + 193e^) [t/S'i sinrei cos(20i — 20o) + ri~^S 2 sin/t2 cos(0i + 02 — 20o)] | . 


(4.8) 


(4.9) 


The substitutions following (^) should be carried out once more to get the secular variation of K 2 - 

The av erage d losse s (^ ), (i.2), ( ^ ) and ( |4.4D (and the similar expression for {dK 2 /dt)), together with the algebraic 
relations ( 2.20 ) and ( |2.2lD provide a complete description of the radiative evolution of the binary system, to order 
in terms of radiative losses. On the other hand the radiative changes of the angles Ki are of relative order 
compared to their secular changes given in ( 2.32|). A n alternative set of evolution equations, in terms of the orbital 
elements a, e and Ki is provided by Eqs. (2.32), ( |T^ ), (^^) and (4.9). It is remarkable that all the contributions from 
the radiative losses of the spins, present in the instantaneous losses average to zero. 


V. CONCLUDING REMARKS 


Our description of the smoothed evolution of a spinning binary system opens up the possibility to study the 
modulation of the gravitational wave forms induced by the eccentricity. Investigations of the smoothed evolution of 
circular orbits have already been presented in several papers ( jl^, and |^). The d etailed analysis of the angular 
evolution equations is a subtle issue. An immediate effect that follows from Eqs. (2.31) of Sec. II is the acceleration 
of the evolution of spin orientations with increasing eccentricity. 

The radiative losses of both the dynamical quantities E and L and of the angular variables Ki and 7 subtended 
by the angular momenta were given here up to order compared to the leading order losses. Among them the 
angular losses of the angles Ki are of order beyond the secular spin-orbit terms given in ( p.32 ). Our previous 
results in I and II are particular cases of th e pres ent radiative loss equations. The one-spin limit arises by ^2 — > 0 
and 01 = 7r/2 (the latter relation stems from ( 2.20| )). For the Lense-Thirring case the additional limit 77 ^ 0 has to be 
taken. We would also like to point out the agreement of the energy and orbital momentum losses with computations 
in a different, noncovariant SSC. There the energy E and orbital momentum L were deri ved from a di fferent action, 
but their radiative losses, given by Rieth and Schafer 0 coincide with those in our Eqs. ( |4.1| ) and ( |4.2| ). However, in 
our approach, the previously unknown radiative evolution of the angles Ki and 7, characterizing the geometry of the 
binary system, could readily been obtained. 
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